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We propose an all-geometric implementation of quantum computation using neutral atoms in 
cavity QED. We show how to perform generic single- and two-qubit gates, the latter by encoding 
a two-atom state onto a single, many-level atom. We compare different strategies to overcome 
limitations due to cavity imperfections. 



The standard paradigm of quantum computation (QC) 
is a dynamical one: in order to manipulate the quan- 
tum state of systems encoding information, local interac- 
tions between low-dimensional subsystems (qubits) are 
switched on and off in such a way to enact a sequence 
of quantum gates. On the other hand, ever since the 
discovery of the Berry's phase, it has been recognized 
that quantum evolutions, besides dynamical contribu- 
tions, can display purely geometrical effects The lat- 
ter, in view of their very geometric-topological nature, 
show an inherent stability against some local perturba- 
tions. It is therefore a natural and intriguing question to 
ask whether one could take advantage of this geometric 
features to the aim of processing quantum information. 
Indeed one would expect the above mentioned robust- 
ness to result in a resilience against some kinds of errors. 
In other terms a geometry-based strategy for quantum 
manipulations is expected to have some built-in fault- 
tolerant features In the context of NMR 1^ and 
Josephson junction based quantum computing g it has 
been show how to use geometrical phases to implement 
a two-qubit gate that, along with the dynamically gener- 
ated one-qubit gates, is universal. 

To achieve an all-geometrical implementation of quan- 
tum computation one is led to consider more sophisti- 
cated, i.e., non-Abelian, structures. This has been orig- 
inally done in Ref. where the so called holonomic 
approach to quantum computation (HQC) has been in- 
troduced. Quite recently a proposal for implementing an 
HQC scheme with trapped ions, feasible with the cur- 
rent technology, has been put forward 0. In this pa- 
per we discuss an implementation proposal for HQC by 
means of neutral atoms in cavity QED. This is to some 
extent related, at least regarding single-qubit operations, 
to the proposal . We shall show how to perform generic 
single-qubit gates by using a single atom. For realizing 
universal two-qubit gate a mapping of a two-qubit state 
onto a single many-level atom will be used. Finally we 
shall propose a strategy to overcome limitations due to 
cavity imperfections. 



I. HOLONOMIC QUANTUM COMPUTATION 

We now briefly recall the basics of HQC d). In the 
HQC paradigm information is encoded in an n-fold de- 
generate eigenspace C of a Hamiltonian H{Xo) belong- 
ing to a d-dimensional parametric family of isodegenerate 
Hamiltonians {H{X)}. The A's represent parameters that 
are supposed to be controllable in the given experimental 
situation. The manipulations of the codewords in C are 
enacted by driving the control parameters along loops x 
in an adiabatic fashion. In this way an initial prepara- 
tion l^o) € C evolves, up to an overall dynamical phase, 
according to the rule |4'o) € C i-^- U{x)a l^^o) where 
U{x)a ~ Pexp J^A is the holonomy associated with 

X by the the M(r7,)-valued connection A = "^-^p- 
One finds § 

Af^irwi^jrw), (1) 

where {|'0"(-^))}S=i denotes an orthonormal basis of 
the degenerate eigenspace C. The set of all possible 
holonomies, obtained by taking all possible x's, is a sub- 
group, known as the holonomy group, of the group U{n) 
of unitary transformations over C. 

When the holonomy group coincides with the whole 
U (n) one can perform universal QC over C by re- 
sorting to geometrical means only. This irreducibility 
condition can be easily stated in terms of the curva- 
ture 2-form F associated with A by the relation F^i, = 
df^Ai, — djyAfj^ — [Af^, Ajj\ jl^. The number of linearly in- 
dependent -F)ji/'s gives a lower bound to the dimension of 
the holonomy group ]l3t . The curvature form F encodes 
for the non-trivial geometric features of the global bun- 
dle of n-dimensional quantum codes over the manifold of 
control parameters. Flat bundles, i.e., with F = 0, have 
no computational power. 

II. DARK STATES IN (N+1)-LEVEL SYSTEM 

Following the previous recipe we will give the expres- 
sion of a parametric Hamiltonian which turns out to be 
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suitable to achieve holonomic quantum computation. Af- 
ter that we will show how such Hamiltonian can be imple- 
mented using neutral atoms in an optical resonator. We 
can, in this way, give a precise physical meaning to the 
abstract objects of the HQC paradigm. The idea rehes 
on the concept of the adiabatic passage via dark states. 
In the simplest case (a.k.a. A-system) we have two states 
(ground states) which are not directly connected, i.e., in 
the Hamiltonian describing the system there are no terms 
which couple these two states - considering an atom in- 
teracting with the electromagnetic field there is no single- 
photon transition between the states we are interested in. 
However they are independently coupled to a third state 
(excited state) in a tunable way. By changing in an adia- 
batic way the coupling constants, it is possible to create 
coherent superpositions of the two ground states and to 
pass from one to the other without populating the excited 
state UM. Our starting point is the Hamiltonian 



H 



N 

E 

fc=i 



^k\e){gk\ + h.c, 



(2) 



which can be seen as a generalization of the A-system 
Hamiltonian and represents a system in which N ground 
states \gk) are coupled to an excited level |e). The level 
structure we have in mind is depicted in Fig. 




FIG. 1. Internal level scheme of the atoms inside the cav- 
ity, with the coupling needed to realize the Hamiltonian Eq. 

(!)• 



By explicitly computing the connection form A and its 
curvature for this system we have checked that it allows 
for universal QC over C for any N . In particular, if 
iV = 2" -f 1 n qubits can be encoded in C. In the next two 
subsections we turn on the cases = 3 and A^ = 5, or 
equivalently n — 1 qubit and n = 2 qubits. We will show 
how to realize singie-qubit rotations - thus, up to a phase 
factor, any single-qubit gate - and a 2-qubit phase-gate 

i. 



A. (3+l)-level system 



One can realize a single-qubit encoding using a system 
described by the Hamiltonian (^) where A^ = 3. As com- 
putational basis we choose the first two ground states, 
i.e., we assign logical values through the identities 



|5i) = |0), 132) = 



(4) 



while the third ground state, jga), plays the role of an 
auxiliary state, which is necessary to achieve every single- 
qubit gate. It is well known that any single-qubit gate can 
be decomposed (up to a phase factor) in the product of 
three rotations, for instance a rotation about the z axis, 
one about the y axis and one again about the z axis, i.e., if 
U is the single-qubit gate we want to build up, there exist 
three numbers, a,/3, 7, such that U = Rz{'-f)Ry{f3)Rz{a), 
the Euler angles. Within our model Hamiltonian a ro- 
tation about the y axis can be obtained by putting the 
relative phases in Eq. (^) (with N = 3) (f)i = 0, and by 
adiabatically changing the amplitudes 0i , 02 ■ In this case 
the connection is just A — A0^dO2 — —isin{0i)ayd02, 
where ay is the y-Pauli matrix. We obtain the unitary 



operator (see App. Al) 



The Hamiltonian admits A^ — 1 dark states \Di), 
i.e., zero-energy eigenstates having no contribution from 
the excited state. The complex couplings flk represent 
our control parameters. Once they are fixed, the cod- 
ing space, i.e., the (A^ — l)-fold degenerate eigenspace 
spanned by the dark states, is the orthocomplement of 
the vector Af-'/' Ek=i^k\gk) {Af = Ek=i W)- % 
setting ^ll — — ■ ■ ■ = ^N-i = 0, we have that the 
coding space C is spanned by the first A^— 1 ground states. 
We write the coupling constants in generalized spherical 
coordinates 



Qi = |rj|sin(6li), 

^2 = |f2|e-''^^ cos(6'i)sin(6'2). 



(3) 



n 



N 



\n\e- 



cosi 



. .C0S[tlN-2 



) sin( 



RyiP) = exp(i/5(Ty), 



(5) 



after a cycle x in the {9i, 6'2)-submanifold, where the an- 
gle (3 is given hy (3 = /g(^) cos{9i)d9id92 and S{x) being 
the surface enclosed by the loop x on the 9i , 6*2 submani- 
fold. Up to a global phase, a rotation about the z axis is 
equivalent to the operator exp(ia|l)(l|), which is easily 
obtained by putting 9i — (j)^ — and by adiabatically 
performing a closed path x' in the submanifold of ^2 , 02 ■ 



The connection is A = 
After a cycle we obtain 



exp(za|l)(l|) =e'"/2^,(a). 



-isin(6l2)a|l)(l|d(/)2 



where a = shi{292)d92d(f>2- 



(6) 
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B. (5+l)-level system 



1. single-qubit gate implementation 



We assume two qubits mapped to a four level system. 
Thus, to implement a two-qubit gate, we need = 5. We 
want to show how to realize a phase gate, which assigns 
a phase only to one out of the four computational basis 
states. As the computational basis we choose \Di) = \gi), 
« < 4, whereby the corresponding coupling constants ^Ik 
are initially set to zero. The logical states can be identi- 
fied as 

l5i> - |00), I52) = |01), b) = |10), \94) = 111), (7) 

the fifth state playing the role of an ancilla. Con- 
sidering a closed path in the two-dimensional sub- 
manifold of the parameter space with coordinates (64, (f)^) 
- the other parameters being kept to zero - the con- 
nection is reduced to the simple form: A — A^jj^dcp^ — 
— i sin^(6'4)|(74)(g4|d(/)5. This gives rise to the holonomy 
(see App. K% 



^Aix") = exp{ia\gi){gi\) 



(8) 



which, according to Eq. (0), precisely represents 
a 2-qubit phase gate, with the phase a given by 
Js(^'') sin{294)d04:d(f>5. 

III. PHYSICAL REALIZATION 

In the remaining part of the paper we discuss a possible 
physical realization of the - up to now quite abstract, but 
also very general - concepts we introduced in the previuos 
sections. Our proposal is based on atoms trapped inside 
an optical resonator (Fig. H). The atoms, which repre- 
sent our qubits, interact individually with laser beams 
and with a single quantized mode of the optical cavity. 





lasers 

FIG. 2. Schematic representation of the setup. 

The manipulation of the qubits (in terms of single- 
and 2-qubit gates) involves only the laser beams, while 
the cavity mode is used to prepare the system every time 
a 2-qubit gate is required. Indeed, as we will see, we 
need to encode the information of two qubits in a single 
many-level atom and this is performed with the aid of 
the cavity mode, to which all atoms are coupled. 



The single-qubit gates are easily implemented if we 
consider (3-|-l)-level atoms, i.e., atoms with three ground 
states coupled by lasers to a single excited state. The 
Hamiltonian of the single atom can be reduced to Eq. (||) 
with N — 3, where the coupling constants flk, fc = 1, 2, 3 
are the Rabi frequencies of the lasers. Therefore for the 
feasibility of single-qubit gates the addressing of single 
levels in single atoms is required. 



2. 2-qubit gate implementation 

Once we have defined qubits, we need a way to couple 
them in a suitable way to make the computation univer- 
sal. In our case we are able to implement any 2-qubit 
transformation in a single atom, having N — b ground 
(or meta-stable) states which can be coupled by tunable 
lasers to a single excited one |l2|. Thus we have to deal 
with a system in which the interacting part is described 
by the Hamiltonian Eq. (0) with N ~ 5 and such that 



we can use the prescription of sec. [IB. By adiabatically 
acting on the coupling lasers, which address the single 
levels of the atom, any 2-qubit operation is achievable. 
The complete picture is based on (5 + l)-level atoms, 
single-qubit information being stored in the single atoms 
and single-qubit operations being performed in any sin- 
gle atom as discussed above, using 3 out of the 5 ground 
states. Let us consider two (5 + l)-level atoms and let 
\a)i and \(3)2 be the (logical) states of the first and the 
second atom, respectively. A 2-qubit gate is performed 
in three steps: 

1. the two-qubit information is stored in the second 
system by the transfer 

|a)i|/?)2-|0)i|O2, (9) 

where a {f3) represents the first (second) digit of ^ 
in binary notation; 

2. since we suppose that any (5 + l)-level atom can 
be driven by the Hamiltonian Eq. (||) , it is possible 
to obtain any (2-qubit) gate by manipulating the 
coupling constants, physically the Rabi frequencies 
ilfe. We have shown above how to obtain the phase 
gate Ua = exp(zQ!|ll)(ll|). In this case it is suffi- 
cient to act only on two of the couplings fJfe, the 
others being turned off; 

3. after the holonomic 2-qubit gate operation, the in- 
verse transformation of Eq. (^ is performed, and 
each qubit is encoded back in one of the two atoms. 

To pursue our purpose what is missing is a method to 
perform the information transfer Eq. (|^). In order to be 
consistent with the holonomic paradigm, the information 
transfer has to be adiabatically performed. 
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A. Information Transfer 

We will suggest two possible approaches to achieve the 
information transfer process. They involve, besides the 
ground states, excited states as well as the single cav- 
ity mode. Thus, such processes will be affected by both 
spontaneous emission from the excited levels and imper- 
fections of the cavity. The first approach, that from now 
on we will call the optical scheme, was envisaged in 
It is based on an adiabatic transfer that leaves almost 
unpopulated the excited levels, thus reducing the influ- 
ence of spontaneous emission. The second one, proposed 
here for the first time, is based on an adiabatic transfer 
that leaves almost unpopulated the cavity mode, thus re- 
ducing the influence of cavity losses. We will call it the 
motional scheme. We turn now to study in detail these 
two different approaches. 

Both of them rely on the level scheme shown in Fig. ||. 
A laser addressing atom i (i = 1,2) with Rabi frequency 
ri'i(i) couples \gi)i to the excited state \ei)i and \g2)i to 
|e2)i, while \g3)i and \g4)i are coupled by the cavity mode 
to, respectively, |ei)i and \e2)i- 
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FIG. 3. Level scheme and coupling laser lights for the in- 
formation transfer. The logical values are also shown (see the 
text). 

Note that the information transfer is equivalent to state- 
swapping between the two atoms, provided the physical 
encoding is the one shown in Fig. ^, i.e., 153)1 = |0)i, 
151)1 = |l)i for the first atom and the 2-bit words, which 
we want to encode in the second atom, are physically 
represented as 1^3)2 = |0)2, 1^4)2 = |1)2, |5i)2 = |2)2, 
[52)2 = 13)2. We have, for the logical state: 

153)1153)2 ^ 153)1153)2 ^ |0)l|0)2 ^ |0)i|0)2 
153)1154)2 -> 153)1154)2 ^ |0)l|l)2 -> |0)i|l)2 
l5l>l|53)2 ^ |53)l|5l>2 ^ |l)l|0)2 -> |0)i|2)2 
l5l)l|54)2 ^ |53)l|52)2 ^ |l)l|l)2 ^ |0)l|3)2 

Note that we used a different idcntiflcation of the logical 
state with respect to sec. 



II B 



To describe the state-transfer process, we will focus on 
the evolution of one out of the two three-level systems 
which are contained in each atom, i.e., the one formed 
by \gi)i, \g3)i and \el)^. 

The evolution of the system, in presence of decoher- 
ence, will be described by a non-Hermitian Hamiltonian 



in the framework of the quantum-jump approach to dis- 
sipative processes 



151 



In the optical approach the following single-atom Hamil- 
tonian is considered: 



H. 



opt 



h[ujeAei)i{ei \ +ujg,\gi)i{gi \ + ujg^\g3)i{g3 \ +ujcb'^b 
(r!',(<)e— *|ei),(5i| + 5fe|ei).(53| + h.c.)] 
*7|ei)(ei| -zK^tfo (10) 



where , LUg-^ and LUg^ are the energies of the state 
|ei)i, |5i)i and \g3)i, respectively, g is the dipole-coupling 
constant between the cavity mode and the atom and b 
is the annihilation operator for the cavity mode. As 
decoherence mechanism we have considered the spon- 
taneous emission 7 of the excited levels \ej)i and the 
cavity loss rate k. We will see that the transfer pro- 
cess is based on a dark state of the compound system 
Hatomi ® Hatom2 <8) Hcavity which docs uot iuvolve any 
excited states and thus, within the adiabatic approxima- 
tion, the main dissipative channel is due to the cavity 
loss rate. 

In the motional approach we assume that the atoms are 
individually trapped in harmonic potentials. The single- 
atom Hamiltonian contains also the harmonic trapping 
potential terms (see for instance [|l8| , ^ ): 

Hl^^ = h[u!eAei)i{ei\ + ujg^\g3)i{g3\ +ujcb^b + va^a 
+ {n\{t)e~'^^'\ei),{g3\+h.c.) 
+ gsin{kx){b\ei)i{g3\ + h.c.)] 
- i-f\ei){ei\-iKb^b (11) 

where a (a^) is the annihilation (creation) operator for 
the harmonic motion, while the sine function describes 
the standing-wave structure of the cavity field, with k the 
wave number of the field and x — [h/{2mi')]^/^{a'' + a). 
We will see that, under certain conditions, it is possible 
to obtain an effective Hamiltonian which involves only 
the cavity mode and the harmonic motion. With such 
Hamiltonian the transfer is based on a dark state with 
respect to the cavity and thus one expects that the most 
important dissipative channel will be the spontaneuos de- 
cay. 

First of all we will show, neglecting any dissipative mech- 
anism, i.e., J — K — 0, that indeed we can obtain the 
state-tranfer by acting in an adiabatic fashion on the pa- 
rameters of the Hamiltonians Eq. (^0|) and Eq. ([Tl]). 
Later on we will carry out in detail the analysis of the 
effects of the dissipative channels. In the next two sub- 
sections we do not consider any decay mechanism. 



1. Optical state transfer 

In the optical approach the starting point is the single- 
atom Hamiltonian Eq. ( p^ ) (where for the moment 
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7, K = 0). Considering the energy hujg-^ as the zero 
of the energy scale and transforming to the reference 
frame described by the Hermitian operator A — ujrb^b + 



WL|ei)(ei| + (WL - U)c)\g3){g3 
gets the Hamiltonian |16[ 



JtA 



in particular the zero-excitation eigen-space is 
spanned by the vacuum state |0)™^|0)2™|0)cav and 
the 1-excitation eigen-space by the three states 

{|i>rio)rio)cav, io)rii)rio)cav, io)rio)rii>cav}, 



H. 



opt 



where is the eigenstate of the free external Hamil- 

^ tonian of the ith atom, it satisfies a|ai|n)f° — ri|n)™. 

fi[-A|ei),(ei| -t- {n ,(t)|ei),(5i| -h56|ei)»(g3| +/i-c.)]-lnside the last sub-space there exists the dark state (dark 

(12) with respect to the cavity mode) 



In writing Eq. (|l^) we have considered the resonance 
condition = {tUg^ — ujg-^) — {ul — Uc) = and we 
introduced the detuning A ^ lul — (o^ei 
the laser frequency ujl and the transition \gi 
The 2- atom Hamiltonian i/opt 
eigenstate 



Tjl 

-"opt 



ujg^ ) between 

) - ki). 

iJgpj admits the 



|*)opt(Xr!'2ff|ffl)l|53)2|0) 

- n\n'2\g3)i\93)2\i) 



^'ig\93)l\9l)2\0)c 



(13) 



Thus by adiabatically changing the Rabi frequencies f2'i 
and applying a "counterintuitive" ( and reference 
therein) pulse (whereby the pulse on atom 2 precedes 
that on atom 1) it is possible to pass from the state 
Iff3)i|5i>2|0)cav to the state |gi)i|53)2|0)cav The state- 
swapping is thus realized. Note that in such a scheme 
(see Eq. (13)) during the information transfer the 1- 
photon cavity state is populated. 



|vl.)ext OC r!'2|l)r|0)r|0)cav ^ ^i' 1 10) 5" 1 1)2" |0)cav • (16) 

In an analogous way to what described in the previous 
Section, by adiabatically changing the Rabi frequencies 
il'i{t) it is possible to pass from the state |l)i™|0)2™|0)cav 
to the state |0)i™|l)2™|0)cav without populating the cav- 
ity mode which, in the case of a nonzero cavity loss rate 
K, is a source of decoherence. 

We are now ready to describe how the state-swapping can 
be performed in three steps starting without photons in 
the cavity and zero-motional-excitation, see Fig. ^ first 
of all the logical state of qubit 1 is swapped onto the mo- 
tional state of atom 1 by adiabatic passage via a dark 
state; then we utilize the dark state ( p^ ) to transfer the 
motional state to the second atom and this is swapped 
onto the internal state (of the second atom) by the same 
adiabatic pulse sequence as in first step, ending with a 
global state without any excitation, i.e., no photons, no 
quanta of harmonic motion. 



2. Motional state transfer 



The single-atom Hamiltonian will be in this case the 
Eq. (11) with 7 = K = 0. Considering the Lamb-Dicke 
limit, i.e., the size of the harmonic trap small compared 
with the optical wave-length, fc[fi,/(2mzy)]^/^ <C 1, and 
writing the Hamiltonian in the reference frame given 
by the Hermitian operator (the zero-energy being ujg^) 
A = ci;i|ei)(ei| -I- ojcb^b + i/a^a one gets 



m 



?i[-A|ei),(ei| + {n',{t)\ei)^{g3\ + h.c.) 



-f 77(a^&|ei)^(53| + abe-'''''\ei)r{93\ + h.c.)] (14) 



where r] = k[h/ {2mv)Y/'^ is the Lamb-Dicke parameter. 
It has been shown in jl^,|l^ that under the conditions 
A 3> \^l'i\^g'q a Hamiltonian which does not involve the 
atomic internal degrees of freedom is obtainable. By adi- 
abatic elimination and using the rotating wave approx- 
imation (RWA) we find the effective Hamiltonian (see 
App. U) 



= h^G^{t)a\b + h.c, 



(15) 



i=l 



where we introduced the coupling parameters 
Gi = —grjVl'iit)/^.. The number of excita- 
tions X)i=i ^I'^i + b'^b is a conserved quantity, 




FIG. 4. Example of information transfer by using exter- 
nal degrees of freedom: \gi,g3) l<?3,5i), photons, within 
the validity of Eq. (^, are never present in the cavity; the 
corresponding logical operation is ll)i|0)2 — » |0)i|2)2. 



IV. DECOHERENCE EFFECTS 

We will now come back to study what are the limita- 
tions imposed by the dissipative channels on the transfer 
process. As stated well before we want to take into ac- 
count excited level spontaneous emission and cavity loss. 
Thus the evolution of the system is described by the Eq. 
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(|10|) and Eq. ( |ll| ) , where the presence of the spontaneous 
emission rate 7 and the cavity loss rate k give rise to a 
non-hermitian evolution. 



A. Optical state transfer 

The optical state transfer is based on the dark state 
Eq. (|l3|). First of all we note that it has a non-zero 
projection onto the 1-photon cavity state and so a lossy 
cavity will tend to destroy such a state. The dark state 
is, indeed, not an eigenstate of the Hamiltonian Eq. ( p^ ) 
for K ^ Q. In order for the influence of the cavity loss 
rate to be small, the condition k Pi-pii{t)dt <C 1 has 
to be satisfied, where Pi-ph{t) — (4'opt|&^&|^opt) is the 
population of the cavity mode during the adiabatic evo- 
lution and T is the time process. Since the integral in 
the inequality is always smaller than Tmaxt{Pi_ph(i)}, 
inserting into the expression Eq. (|l^) Gaussian-shaped 
laser pulses with peak value f2, time separation aT and 
variance one can give the following requirement 



KTmax({Pi_ph(i)} = kT [l + 



< 1. 



(17) 



Furthermore, in any real process (finite time, finite en- 
ergy), the state of the system will precess around the 
dark state, instead of following it in a perfect adiabatic 
way. This means that some population reaches the leaky 
states \e)i. One should calculate the population of 
such unwanted states during the adiabatic process and 
impose the condition 7 Pe{t)dt <C 1. We give here 
and in what follows some simplified conditions. The con- 
dition for adiabatic passage can be stated as 



r2r> 1. 



(18) 



Here Hfl is a parameter expressing an estimate of the 
global value of the differnce between the dark-state 
energy (zero) and the smallest (non-zero) eigenenergy, 
Emin{t), of the Hamiltonian Eq. (|l^ ) during the transfer 
process. The condition on 7 reads; 



< 1. 



(19) 



The transfer time, T, is bounded from below and from 
above. Let us compare the on-resonance case, i.e., A = 



and the far off-resonance case, i.e., 
evaluate as a time average of 
regimes we find 



A >_[£]', I, g. We 



(t) |11|. In these 



1 



"cff 



1/2 



with fl^g 



Since we know that the adiabatic passage works better 
when a ~ T [|n|, from Eq. (|l^, we learn that it is 
favourable to have g ^ il. The tranfer time is restricted 
to be: 



— ^ > r > — - — 



C, 



(22) 



where C = 1 for A = and C = 8A^/{fl'^ff) for large de- 
tuning. Provided that the previous inequality is satisfied 
on both sides by a factor a, we obtained that, at least, it 
must be K7 < {g/aY for A = and K7 < {g / a)^ {fl / A)'^ 
for large A. Then it is clear that the optical scheme works 
better in the on-resonance regime. 



B. Motional state transfer 

For the transfer procedure via motional state swap, 
it is possible, as we have seen, to obtain a Hamiltonian 
involving just the external degrees of freedom and the 
cavity mode. This is true also in the presence of deco- 
herence processes, where one obtains the non-Hermitian 
Hamiltonian 



A- 



iKb^b - 



47 



'A2 + 72 



(23) 



The details of the calculation are given in App. B S , 
where also the Hamiltonian without RWA is shown. To 
obtain Eq. ( p3[ ) the conditions A ^ \fl'i\,rig are im- 
posed. If T is the process time, the conditions we can 
give on 7 are (first order pertubation theory): 



7T 



< 1. 



(24) 



The motional scheme is based on a dark-state with re- 
spect to the cavity mode (see Eq. (p^)). The actual state 
of the system during the evolution will precess around the 
dark state and so we will have a certain population Pi-ph 
of the 1-photon cavity state. In order to give a condition 
on K we have to estimate such population along the adi- 
abatic process. For Gaussian-shaped pulse we found (see 
App. 0): 



^'i-ph < — e 

T 



V A 



T 



(25) 



The condition on k reads: P/iPi-ph <C 1. Thus consid- 
ering for convenience m = maxjr/g, fi}, we see that the 
^20)process time has to be (a ~ r): 



(21) ^ f^V>P> ^ ^ 

7 \mj {tviY g2 



(26) 
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First of all note the inverted role of 7 and k in the opti- 
cal and in the motional scheme in restricting the process 
time. Then, if the inequality is satisfied on both side 
by a factor a we have that n'y < {g / of {tt^Q, / . We 
see that, in the motional scheme, we have a much more 
strict restriction on K7 due to the presence of the small 
Lamb-Dicke parameter. 



C. Modified optical sciieme 

In this section we describe how to enrich the optical 
scheme to obtain a new scheme, where the role of the 
decoherence mechanisms is the same as in the motional 
scheme, without the Lamb-Dicke parameter. 

In the case of large detuning , i.e., A » |ri'i|, g, we can 
perform an adiabatic elimination of the excited states. 
Starting from the Hamiltonian Eq. (p^, we obtain the 
following 2-level Hamiltonian (sec App. BJ^ and fl^): 



(30) 



The condition on k is TnPi^ph ^ 1. The presence of the 
auxiliary levels impose on 7 the restriction 



< 1. 



(31) 



Thus, for the process time, T, we have (as usual a ~ r) 

1 /A\' 1 K /AV , , 

-(tt) . (32) 

which is the same condition we gave for the motional 
scheme, but without the small Lamb-Dicke parameter. 
For instance the condition on 7K, obtained as explained 
before, reads 7K < {gjaf'T^ . 

D. Summary 



en ,1 



+ 



A -I- 17 
A + 17 
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A -I- 17 



(&^l53)»(gi| + 6|5i).(53|)-*^fot6. (27) 



The total Hamiltonian 



TjOpt 



^X.i^ when 



K = 0, admits, as expected, the zero-energy eigenstate 
Eq. (p^. If it is possible to compensate the total Stark 
shift of the ground state \g\)i - i.e., the real part of the 
complex energy of \g\)i in Eq. ( p7| ) - one obtains a Hamil- 
tonian which has a similar structure to the one obtained 
for the motional scheme, Eq. (|l5|): 



A2 -I- 72 



A -I- 17 



g^\ 
A + 17 



(&^l53)»(gi| + 6|5i).(53|)-*«fotfe. (28) 



Physically the compensation can be realized by coupling 
the state j^i) with another auxiliary state via a laser 
pulse which has the same Rabi frequency as the one that 
couples \g\)i to |ei)i, i.e., Sl'^, and a detuning —A. For 
simplicity, the auxiliary state is supposed to have the 
same spontaneous emission rate 7 as have. When 
7 = (while n, can be non-zero), the Hamiltonian Eq. 
( p8| ) admits the dark state with respect to the cavity 
mode 



I*) CX r!^|5l)l|53)2|0)cav + f^'l|<?3)l 151)2 |0)c 



(29) 



To obtain the requirement on k, we estimate the popu- 
lation of the 1-photon cavity state during the adiabatic 
transfer. Using the same notation we used above we ob- 
tained (sec Appendix 1^): 



In summary, using the motional scheme it is possi- 
ble to reduce the effects of k, by increasing the process 
time, since it is based on the adiabatic transfer via a 
dark state with respect to the cavity. The process can- 
not be too long in order to avoid spontaneous emission 
effects, which are however reduced by choosing a large 
detuning. On the other hand the optical scheme, based 
on a dark state with respect to the excited levels, can 
avoid the effect of the spontaneous emission by increas- 
ing the process time. In this case, however, the process 
must be fast with respect to the inverse of the cavity loss 
rate. Finally, provided we can use a stabilizing pulse to 
compensate for the time-dependent energy of the ground 
state (see also ) , the optical scheme in the large detun- 
ing regime could be able to operate in the same fashion 
as the motional scheme but faster (basically by a factor 
of the order of the Lamb-Dicke parameter 77) and then to 
allow for a larger spontaneuos emission rate 7. 



E. Numerical results 

We performed numerical simulations for the suggested 
methods. They confirm what we have stated in the last 
sections. We show in Fig. (^) and Fig. ^ the transfer 
fidelity, respectively, for the optical scheme and for the 
motional one, under the same condition of the parameter 
7, K, f2 and g. We expressed all the quantities in unity of 
the coupling, g, between the atoms and the cavity mode. 
We simulated the evolution of the system as driven by 
Gaussian-shaped pulses. We have followed the evoltion 
of the projections of the state onto the bare states, i.e., 
\gi) and je^). We are interested in the population of the 
transfer's target state Fig. (||). 
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Optical scheme - T = 6 x lO^/.g and k = lO-^g-.-lO- 





time (los) 

FIG. 5. Adiabatic evolution of the system. The figure rep- 
resents the evolution of the projection of the dark-state, for 
the optical scheme, onto the initial basis states. We see that, 
essentially, only the projection onto the states involved in the 
dark-state are populated. 



It is this population that it is plotted, as a function of 
AC and 7, in the Figs. and (|^), the time is fixed. In 
such shown figures we use: 57 — 0.05(7, A = 0, for the 
optical scheme, A = lOg and 77 = 0.1 for the motional 
scheme. 




Optical scheme - T = 6 x W/g and y = 10-^g...lO-^g 



03 
J3 



O 




FIG. 6. Population of the "target" state for the optical 
scheme as a function of the spontaneous decay rate, 7, for 
various values of the cavity loss rate, k (up); as a function of 
the cavity loss rate, k, for various values of the spontaneous 
decay rate, 7 (down). 



As expected from the previously given inequalities the 
optical transfer is more affected by the presence of the 
cavity loss, than by the presence of the spontaneous emis- 
sion, the vice versa is true for the motional scheme. In 
this parameter regime, the optical scheme works much 
better by means of both fidelity and transfer time, than 
the motional one. However is worth to remind the differ- 
ent role played by the two decoherence channels as it has 
been summarized above. In principle for a bad cavity and 
a small spontaneous emission rate, one shall have, wait- 
ing enough time, a fidelity unreacheable with the optical 
scheme. 
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Motional scheme - T 



a, 
o 
a, 




Motional scheme - T = 6 x l(f / g and k = 10 ^5. ..10 




FIG. 7. Population of the "target" state for the motional 
scheme as a function of the spontaneous decay rate, 7, for 
various values of the cavity loss rate, k (up); as a function of 
the cavity loss rate, k, for various values of the spontaneous 
decay rate, 7 (down). 



V. CONCLUSION 

In this paper we have presented an implementation of 
quantum computation based on holonomic operations, 
i.e., non-Abehan geometric phases, with trapped neutral 
atoms. The method is a generalization of a previous pro- 
posal exploiting ions |^]. We showed explicitly how to 
realize a general unitary transformation on n qubits en- 
coded in a (2" -I- 2)-level atom. We also developed a 
scalable scheme, where each qubit is stored in one atom. 
In this case, two-qubit operations are performed in one of 
the two atoms, which encodes both atoms' initial logical 
state after a properly designed transfer process. We dis- 
cussed in detail two possible procedures for such a state 
transfer, relying on adiabatic passage via a dark state of 



the two-atoms-plus-cavity system. The two schemes (the 
second of whose is originally proposed here) differ in the 
atomic degrees of freedom involved - respectively, inter- 
nal or external. We discussed advantages and limitations 
of both schemes in the presence of decoherence, for dif- 
ferent parameter regimes, and we found that the second 
proposal is more suitable for a situation with a compar- 
atively lossy cavity, as it might be the case, e.g., with 
atomic micro-traps coupled to surface-mounted micro- 
cavities in the context of the so-called Atom Chips 20|. 
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APPENDIX A: HOLONOMIES 
1. Single-qubit gate 

As described in the text any single-qubit gate can be realized by holonomic means if we consider the Hamiltonian 
Eq. (^) with = 3. We write the coupling parameters using a kind of spherical coordinates: 

fii = |f]|sin(6'i), 

^2 = |J]|e-''^= 008(6*1) sin(6'2), 

^3 = \n\e-'''"> 008(6*1) 008(6*2). (Al) 

The Hamiltonian admits two zero-eigenvalue eigenvectors jV'^) and \ip'^), which can be written in terms of the ground- 
states |(7fe), A; = 1,2,3 as 

IV') = cos(0i)|5i) - sm{9i)[e-'*' sin(02)|g2) + e"^*^ cos(02)|ff3)] 

1^2^ = _e-^4'2 cos(02)|ff2) + e-'^' sin(02)|53). (A2) 

We fixed the two relative phases 02 and 03 to zero and calculated the connection components = {ip^ld/dX^lip^) 
where {A^} = {9k, (pj'jk = 1, 2, 3; j = 2, 3}. We obtained the connection 

A = Ae^de2 = -isa\{ei)ayde2, (A3) 

where Uy is the y-Pauli matrix. The related unitary operation is U{x)a — Pexp ^_^A — Pexp Ag^dd2, where the 
integral is along a loop x in tti6 sub-manifold 9 1,62- The line integral can be converted using the Stokes theorem in 
the surface integral /g^^) 62*^^1 "^^2, where 

Fe^e2 = de^Ag^ - dg^Ag^ = -i cos(6'i)crj, (A4) 

are the components of the so-called curvature 2-form, and S{x) is the surface enclosed by the loop x in the {61, 6*2)- 
plane. The unitary operator takes the form: 

[/(x) = exp-icTy / cos{0i)d9id92 (A5) 

with cFy the y-Pauli matrix, i.e., we have obtained a qubit rotation around the y-axis of the Bloch sphere. In the same 
way one can obtain the gate Eq. (^, fixing 0^ = cf)^ — 0. In this case \ip^) = j^i) for any values of the parameters and 
so, after a cycle in the remaing parameter submanifold, only the state 152) acquires a phase. We stated also in the 
text that by choosing in a suitable way the loop in the manifold {0fc,0j;fc = l,2,3;j = 2,3} it is possible to obtain 
any single-qubit gate. Showing the feasibility of the 2-qubit phase gate is enough to conclude the universality of this 
approach. We stress once more that in this approach, anyway, instead of thinking how a particular gate is decomposed 
in single- and two-qubit gates and then to try to realize them, can be much easier from the experimental point of 
view to search suitable loops (experimentally it means to search for the simplest loops) to perform the quantum gate 
we want. Indeed we have that, given any unitary operator, U, there exists a closed path, in the parameter space such 
that the holomony it generates coincide with U. 
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2. Two-qubit gate 



Realizing the phase-gate does not require more effort than realizing a single-qubit rotation. In this case the 
Hamiltonian is given by Eq. (^) with = 5. We write the coupling parameters using a kind of spherical coordinates: 

Qi = |r2|sin(6ii), 

^2 = \n\e~"''^ cos(6'i)sin(6l2), 

^3 = \n\e~''f'^ cos(6'i)cos(6'2)sin(6'2), 

^4, = \n\e^''f'^ cos(6'i) cos(6'2) cos(6'3) sin(6'4), 

= \n\e~''f"- cos(6'i) cos(6'2) cos(6'3) cos(6'4). (A6) 

This Hamiltonian admits four zero-eigenvalue eigenvectors \ip'') {i = 1,2,3,4), which can be written in terms of the 
ground states l^fe) (fc — 1, 2, 3, 4, 5) as 

IV') = cos(0i)|<7i) - sin(0i){e-''^= sin(02)|ff2) + cos(02) [e-"^'\93) + cos(f?3) (e"^"^^ sin(04)|54) + e^^'^^ cos(04)|.95))] }, 

= e-'^' cos(02)|52) - sin(02) [e"^*^ sin(03)|53) + cos(03)(e-'^^ sin(04)|54) + e~'^'' cos(04)|55)] , 
1^3) = e-'*^ cos(03)l53) - sin(03)(e-'** sin(04)|<?4) + e"^^^ cos(04)|<?5)), 

1^4) ^ _g-.04 cos(04)|54) + e-'*^ sin(04)|<?5). (A7) 

It is worth noting that, when all the parameters (actually the angles 0i) are fixed to zero, the previous eigenstates 
coincide with the 4 ground states .., 1^4): we will consider always paths which start and end at such a point. We 
fixed to zero the parameters 01,^2,^3 ~ i-e., the coupling constants ^11,^2,^3 ~ and the relative phases 02, 03, '^4, 
thus the connection gets the simple expression 

A = A^J(j)5 = -isin2 6i4|54)(.g4M(/>5- (A8) 

It gives rise to the unitary operator U{x")a — Pexp f Ag^d62, which once again can be easily calculated by using 
the Stokes theorem to convert to a surface integral the'line integral. Indeed, by introducing the curvature 2-form F, 
which in this case has just one non-zero component F^^g^ = isin(204), one gets 

[/(x") = exp-z / sin(204)d04d05|54)(54|. (A9) 
Jsix") 

The last expression is precisely a phase gate, i.e., an operation that assigns a phase - equal to the surface-integral - 
to one (134)) out of four states. 



APPENDIX B: EFFECTIVE HAMILTONIANS 

1. Three-level atom in a cavity: Effective Hamiltonian 

In this section we are going to find out an effective, approximate Hamiltonian, starting from a well-known Hamil- 
tonian in quantum optics. We consider a single three-level atom interacting with a cavity and we want to take into 
account also the dissipative terms, i.e., as in the previous Section, the spontaneous emission from the excited level 7 
and the cavity loss k. Such information can be embodied in the effective Hamiltonian ( |To| ) as follows: 

i/opt = n[-{A + Z7)|e)(e| - t^b^ + {n{t)\e){gi\ + gb\e){g3\ + h.c.)]. (Bl) 

We put the state vector \tp)£C^ Hcav of the system in the form 

lil) = ^(Ce,m|e,TO) + Cg^.jn\gi,m) + Cg^^rn\93, m + 1)) (B2) 
m 

and by the Schrodinger equation we found the equation of motion for the coefficients Cg-^^m, Cg^^m and Ce^m- 
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'i-Cgs,m = gVm + lce,m " «K(m + l)cg3,„. (B3) 
From the first equation one gets: 

Jo 

then substituting the coefficient Ce,m in the two fast expressions of the Eq. ( |B3| ) with the previous expression, imposing 
that Q,g\/m + 1 <C ^^A^ + 7^ ^nd neglecting the terms of higher order in fi/^/A^ + 7^ a,nd g^/rrT+T/ y^A^+t^ 
one eventually obtains the equations: 



A + Z7 ^ 7 A + 17 

These equations can be equivalently derived starting from a 2-level system - with internal states \gi) and \g^) 
interacting with the cavity mode by the effective, approximate Hamiltonian (]27|). 



2. Trapped atom in a cavity: Effective external Hamiltonian 

Let us write the non-Hermitian effective Hamiltonian, i.e., considering also the dissipative terms, for a harmonically 
trapped 2-level atom inside a cavity QED: 

H={oJo~ i-f)\e) (e| + (wc - iK)b^b + t/flta + (f7|e) (5|e*"^* + h.c.) + r;(at + a){gb\e) {g\ + h.c). (B6) 

Here, cjq is the energy difference between the ground and the excited atomic level, Uc and h are respectively the 
energy and the annihilation operator for the cavity mode, v and a are the energy and the annihilation operator for 
the harmonic motion, il and ujl are the Rabi frequency and the frequency of the laser light, 77 is the Lamb-Dicke 
parameter, g is the dipole cavity-atom coupling constant and 7 and k are respectively the decay rate from the excited 
state and the cavity loss. 

From now on we do not consider in the equation the cavity loss, because its effect on the final Hamiltonian is trivial 
- indeed at the end it adds the term —iub'^b: so we re-introduce it just in the final Hamiltonian(s). Making the 
canonical transformation a ae^*''*, b — + be"^'^"* and writing the state vector jV') G (g) Tioxt ® Ticav of the system 
in the following form: 

1-0) = '^{ce,n,va\e, n,m-l) + Cg^n,m\g, n, to)), (B7) 

n,m 

we found the equation of motion for the coefficients Ce,n,m and Cg^n,m {g, ^ real): 



«Ce,„,m = (t^o - «7)ce,n,m + f^e *"-^*Cg,„,„i_i -|- jygc (^/r^me*''*Cg,„_l,m + V(n + l)me ''"^'Cg^n+i^-m), 

iCg,n,m = Oe^^-^'ce.n^m+l -|- ?7ge"^"* ( V^^^e'^'ce,™-!,™ + -f l)me"*''*Ce,„-|-l,m)- (B8) 

Integrating the first equation yields 

t 

e~ 







-(-o-.7)(*-r)[j7e— .rcg,„,™„i + 775e-'"="(v^e'^^Cg,„_i,™ + ^{n + l)me-^^^c<,,„+i,„)]dr (B9) 



then integrating by parts, considering the resonance condition ujl + i^ — and introducing the parameter A = w^ — 
one gets the expression 



e 



(Slcg, + r]g{y/nmcg^n-i.m + \/(ri + l)rne ^*'''^Cg,„+i,m)) 







e HA+J7)r _|_ rig{^/mncg^n-i,m + \J [n + l)me ^"'^c<,,„+i,m))dT 



A + ij Jq dr 
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If in the previous expression we keep the terms up to the first order in fJ/ ^ + 7^, r\gl \/ + 7^, we obtain 



3 — lCj£,t 



-(f^Cg,n,m-l + VaiVnmCg^n^i^ra + V {n + l)me ^'"'^Cg^n+i^m)). (BIO) 



(A + *7) ^"'^^^^ 

We can now substitute such an expression in the equation for Cg^n.m and to first order in 77 we find 



A + n 

r]gQ I — r- — 

+ ■ , . {\jn{ra + l)cg,„_i,m+i + V ('^ + l)nT-Cg,„+i,m-i) 
A + Z7 

+ '^-^^ (e-^^'^V(" + + l)cg,„+i.„+i + e^"^*v^Cg,„-i,^-i). 
A + 47 



(Bll) 



The first term contains a Stark shift, namely Ari^/(A^ +7^), and a decoherence part 7ri^/(A^ + 7^); the other terms, 
while leaving unchanged the internal atomic states, involve the external (harmonic) atomic states and the cavity state. 
Thus it is possible to write down an effective Hamiltonian for the external atomic degrees of freedom and the cavity 
QED: 

U,^ = -^l^iab^ + e2-*at&t + h.c.) - inh^h ~ i^f^. (B12) 
A + 17 A^ + 7^ 

If the RWA is applicable (i.e., in this case v 1/t where t is the time of the process we are interested in) the previous 
expression takes the form Eq. (|2F 



RWA ^ ^(«6t + „t;,) „ ,,6t6 „ (B13) 
A + 17 A^ + 7^ 



APPENDIX C: REMARKS ON THE ADIABATIC APPROXIMATION 

The optical scheme after having introduced the pulse to compensate the light shift of the ground state \gi) and the 
motional scheme rely essentially on the same Hamiltonian. Thus we will study here the generic problem of adiabatic 
transfer in a system driven by such a parametric Hamiltonian and we will eventually consider the actual expression 
of the parameters involved in the specific scheme. The starting point is the 3-level coupling Hamiltonian 



H = D\i){2,\ + ;i(Gi(t)|l)(3| + G2(t)|2)(3| + h.c), 



(CI) 



where we have explicitely shown the time dependence of the coupling constants Gi, i — 1,2. Note that for the motional 
scheme D = 0. The eigenfrequencies of the Hamiltonian Eq. (CI) are: Eq ~ 0, E± ^ D ± \/D^ 



G^ff, where we 



defined GcS — \/ G\. The respective eigenstates can be written in the foUwing form: 

\i,^{t))^GZ^{G^\\) 



IV'±(t)) 



{El 



-Gi|2)), 
)-i/^(Gi|l) 



-G2|2)+i?±|3)). 



(C2) 



Troughout the paper we have called the zero-energy eigenstate jV'o) the dark state. We prepare the system in the 
state |1). If at the beginning of the process G1/G2 = 0, we have that the dark state concides with such state. Then 
we suppose that the coupling constant Gi are modified slowly - adiabatically - towards the ratio G1/G2 00. The 
adiabatic theorem |2^ tells us that the state of the system preceeds around the istantaneous eigenstate |?/'o(i)), the 
asimptotic state being the state |2). The population of the other (two) states, P±, can be in principle evaluated in 
the adiabatic approximation byusing the expression 



P±{t) = 



iI>±{t) 



— ?Ao(t) )dT 



(C3) 



where a±fi{t) — j^{E±{t') ~ Eo{t'))dt'. An upper limit of the population is |^ 
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p± < 



{i:±\dH/dt\^o) 



Using the various expressions described before, one obtains 



{il^±\dH/dt\4,a) 
h^{E± - Eor 



= 1 



'-'Is 



1± 



D 



-3/2 



G1G2 — G1G2 



(C4) 



(C5) 



Considering _D = and Gaussian-shaped couphng constant Gi = Gexp[— (t/T — of /t^] and G2 — G'exp[— (t/T 
aY /tY in the previous expression we obtained for both states 



P+ < P ~ max < — „ „„ . sech — ^ 



(C6) 



The r.h.s. of the previous inequahty diverges at t — > ±00. We consider only a finite interval - the same we used in 
the numerical simulation of the transfer process ~ where anyway the condition on the initial and final value of the 
ratio G1/G2 are pretty well satisfied. In such an interval the r.h.s. takes its maximum value when Gi = G2, i.e., at 
the center of the interval (see also the numerical simulation). One obtains 



P+ < 



(C7) 



When D ^ the population for the two unwanted states is different. 

In the text we are interested in the population of the leaky states to give a condition on the rates 7 and k. In the 
generic system considered here, this amounts (as we will see soon explicitly) to evaluate the population P3 of the state 
1 3). From the expression for the eigenstate of the Hamiltonian we see that 



E^ 



Gl 



El 



= (V'+I^) 



E. 



Gl 



El 



(C8) 



where is the state of the system during the adiabatic evolution. For the optical scheme we have the identification 
|1) = l5i)i|53)2|0)cav, |2) = 133)1 l5i>2|0)cav and |3) = |53)i|53>2|l)cav for the states and G = gl^/A for the coupling 
constan t pe ak value. In this case D = g^/A 7^ 0; we have estimated that, in our simulation, the term containing D 
in Eq. ( |C5|) is of the order of the unity for P_ and much smaller than 1 for . Th us a n estimate of P3 can be given 
as in the case D — Q. In such a way, and substituing the expression for G in Eq. (C7), we obtained Eq. (|3^). For 
the motional scheme we have the identification |1) = |l)i|0)2|0)cav, |2) = |0)i|l)2|0)cav and |3) — |0)i|0)2|l)cav for the 
states and G — grjil/A for the coupling constant peak value. Note that in this case D = 0. Substituing G in Eq. 
( |C7|) we obtained Eq. (||). 
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